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ABSTRACT 

This  paper  presents  a new  method  for  assessing  the  residual  strength  of  structural  components  which  have 
sustained  impact  penetration  damage.  The  method  is  based  upon  the  application  of  the  Weibull  probability 
density  function  to  account  for  the  significant  scatter  in  residual  strengths  which  result  from  the  random 
nature  of  impact  damage  detail  (cracks,  holes,  and  spall  surfaces).  To  measure  the  extent  of  an  impact 
damage  pattern,  a new  parameter  is  presented  which  s based  upon  the  solution  for  an  inclined  crack  in  a 
uniaxial  tension  field  and  an  experimentally  obtained  fracture  criterion  for  mixed  mode  crack  configurations. 
Test  results  are  presented  for  residual  strengths  of  machined  damage  specimens  to  ascertain  the  statistical 
nature  which  arises  from  extreme  damage  detail,  smooth  holes  and  sharp  cracks.  Finally,  probability  of  failure 
curves  are  given  for  both  machined  (7075-T6)  and  impact  damage  (7075-T6and  2024-T8l)in  aluminum  uniaxial 
tension  panels. 


I.  Introduction 

The  possibility  of  impact  damage  during  the  operational  life  of  a structure  can 
establish  the  necessity  for  a fracture  mechanics,  residual  strength  assessment  even 
when  accurate  NDT  and  inspection  procedures  indicate  highly  safe  performance 
under  ordinary  operating  conditions.  Such  damage  can  obviously  arise  for  structures 
operating  under  the  risk  of  impact  from  natural  or  ballistic  projectiles.  Other  si- 
tuations exist  where  impact  damage  can  occur  from  common  operational  accidents 
such  as  a back-hoe  striking  a pressurized  pipe  line  while  digging  in  an  open  field. 
Various  examples  throughout  the  literature  [1-6]  indicate  that  design  engineers  are 
aware  of  the  problem  and  have  attempted  to  include  rational  consideration  of  impact 
damage  in  the  design  process.  The  past  approaches  as  summarized  in  [7]  have  been 
attempts  to  estimate  the  residual  strength  of  an  impact  damaged  structure  from  a 
deterministic  mechanics/empirical  viewpoint  by  utilizing  simplified  concepts  from 
conventional  fracture  mechanics.  Thus,  procedures  have  evolved  to  assess  the  degree 
of  damage  in  terms  of  some  “effective  flaw  size"  and  relate  it  to  a critical  operational  load 
through  an  “impact  damage  fracture  toughness”. 

This  paper  presents  an  approach  which  acknowledges  that  impact  damage  is  an 
extremely  complex  degradation  of  a structural  component.  For  a given  set  of  impact 
conditions:  projectile,  velocity,  angle  of  impact,  structural  configuration,  and  ma- 
terials, the  detailed  nature  and  extent  of  damage  is  beyond  prediction  by  the  current 
state  of  solid  mechanics.  Moreover,  because  of  the  apparently  random  nature  of  the 
damage  detail  (cracks,  hole  ;,  and  spall  surfaces),  the  proposed  approach  for  the 
evaluation  of  residual  strength  is  based  upon  a statistical  rather  than  a deterministic 
foundation.  In  this  manner  a statistical  probability  distribution  is  utilized  to  account 
for  the  scatter  which  is  apparent  in  the  residual  strengths  of  impact  damaged 
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components  as  a result  of  the  insensitivity  to  the  damage  detail  of  available  damage 
measures.  A methodology  is  established  wherein  one  employs  a new  fracture  me- 
chanics based  parameter,  dim , to  estimate  the  extent  of  damage  from  a specific  impact 
penetration;  then  through  a corresponding  statistical  function  determines  the  prob- 
ability of  failure  as  a function  of  the  operational  load  for  a structural  component 
containing  dim. 

The  following  sections  present  the  basis  for  d,m,  application  of  the  Weibull 
distribution  function  to  correlate  the  probability  of  failure  with  d,m.  and  a demon- 
stration of  the  probability  approach  in  conjunction  with  residual  strength  data  for  high 
strength  aluminum  specimens  containing  machined,  controlled  damage  (cracks  and 
holes)  and  ones  containing  ballistic  impact  damage. 

2.  Impact  damage  measure 

Before  an  evaluation  of  residual  strength  can  be  made,  it  is  necessary  to  measure  the 
extent  and  nature  of  damage  resulting  from  a given  impact  penetration.  Figure  I 
illustrates  a typical  impact  damage  configuration  as  might  result  from  an  impact 
penetration.  Note  that  in  general  a complex  damage  detail  arises  from  which  a 
resulting  mechanism  of  fracture  propagation  will  initiate  when  a critical  load  is 


Principal  Load  Direction 


Hole. 


Spall  surface 


Perimeter  of  Impact 
Oamaqe  Zone 


djm  • max  (p  sin'  $ ) 
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rigure  I.  Measures  of  impact  penetration  damage 
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reached  in  the  structure.  References  [8  and  9]  provide  definitions  of  the  two  hereto- 
fore established  measures  of  impact  damage.  The  first  is  the  lateral  damage,  LD, 
defined  as  the  greatest  dimension  characteristic  of  the  extent  of  detectable  impact 
damage  including  through-hole,  cracks,  and  spall  surfaces.  This  measurement  is 
independent  of  load  orientation.  In  contrast  the  second  measure,  termed  the  trans- 
verse lateral  damage,  TLD,  is  the  maximum  damage  as  measured  perpendicular  to  the 
direction  of  principal  loading. 

TLD  appears  to  have  been  the  more  successful  parameter  for  correlation  with 
residual  strengths  from  all  reports  in  the  literature.  However,  one  can  observe  that 
TLD  as  well  as  any  simplified  damage  parameter  is  insensitive  to  the  smoothness 
(hole)  versus  the  sharpness  (cracks)  which  may  result  from  a given  impact  damage 
situation.  Even  for  the  cases  where  cracks  are  certain  to  be  associated  with  impact 
damage.  Fig.  2 illustrates  a weakness  in  TLD.  Here  is  shown  a series  of  extreme  cases 
of  damage  detail  where  cracks  are  oriented  at  varying  angles  to  the  principal  load;  the 
critical  load  is  shown  to  differ  greatly  with  crack  angle  in  test  data  described  in  a later 
section;  yet  TLD  is  constant.  To  enable  some  influence  of  this  orientation  of  damage 
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to  be  accurately  incorporated  into  an  impact  damage  measure,  the  following  new 
measure  has  been  devised. 

The  stress  intensity  factors  for  a crack  inclined  at  some  angle,  p,  (Fig.  2b)  to  the 
principal  load  direction  in  uniaxial  tension  is  known  as  [10] 

K|  = o-Vrrp/2  sin2  /3  ... 

K„  = crVirp/2  sin  p cos  p 

where  p is  the  total  crack  length.  Fracture  tests  of  structural  materials  under  mixed 
mode  conditions  [11,  12]  show  that  at  least  for  angles  of  P ranging  from  twenty  to 
ninety  degrees,  the  following  relationship  provides  a reasonable  fracture  criterion  for 
many  applications. 

Kf+Ktf-Cj  (2) 

where  C„  is  the  conventional  fracture  toughness. 

Substitution  of  (1)  into  (2)  leads  to  the  following  relationship  between  fracture 
stress,  ov,  and  damage  parameters  p and  p. 

ac  = C0V2Mp  sin2  pY'  (3) 

This  equation,  which  was  derived  for  a single  inclined  crack,  suggests  a parameter, 
dim , which  can  be  extended  for  use  as  a measure  of  impact  penetration  damage. 
Specifically,  dim  is  defined  to  be  the  maximum  value  of  the  term  in  parentheses  in  Eqn.  (3) 
for  the  set  of  all  possible  inclined  cracks  which  can  be  constructed  within  the  impact 
damage  zone. 

dim  = Max  (p  sin2  p)  (4) 

As  illustrated  in  Fig.  I,  to  measure  dim  from  a given  impact  penetration  damage 
pattern,  a perimeter  is  first  established  around  the  detectable  zone  of  damage:  cracks, 
hole  and/or  spall  surface.  The  impact  damage  measure  becomes  the  maximum  value 
of  a damage  chord  times  the  square  of  the  sine  of  the  angle  between  the  chord  and  the 
principal  load  direction  as  all  chords  characteristic  of  the  damage  pattern  are  ex- 
amined within  the  perimeter. 

Note  that  for  most  instances  of  normal,  as  opposed  to  oblique,  impact  damage  in 
nearly  isotropic  materials,  the  resulting  damage  is  fairly  radially  symmetric  and  often 

TLD  = Max  (p  sin  p)  = dim  (5) 

where  p ^ 90°.  For  that  case  a good  statistical  correlation  between  residual  strengths 
and  TLD  can  be  anticipated. 

However,  for  cases  of  highly  skewed  damage  patterns  resulting  from  an  oblique 
impact  or  a highly  anistropic  material/damage  response.  dim  more  accurately  measures 
the  extent  of  damage  for  use  in  a residual  strength  determination.  It  is  postulated  that 
while  neglecting  some  possible  crack  damage  interactions  not  accounted  for  in  Eqn. 
(D.  dim  includes  a rational  variation  with  damage  orientation  and  relates  to  the 
fracture  stress,  <rf,  as  a square  root  “flaw  size”  inherent  to  conventional  fracture 
mechanics  toughness  parameters.  This  makes  it  an  acceptable  parameter  for  the 
statistical  correlation  of  residual  strengths  of  impact  penetration  damaged  com- 
ponents, especially  when  the  strengths  are  crack-controlled. 

Note  that  dim  does  not  include  parameters  which  reflect  material,  plate  thickness 
or  certain  aspects  of  damage  detail  such  as  sharpness,  or  hole  versus  crack  size.  All  of 
these  factors  influence  residual  strength.  The  influence  of  structural  parameters  such 
as  material  and  thickness  can  be  isolated  by  selective  grouping  of  residual  strength 
data  before  statistical  analysis  is  performed.  The  influence  of  damage  detail  not 
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included  in  dlm  is  accounted  for  in  the  spread  of  the  statistical  distribution  of  residual 
strengths. 


j 


3.  Residual  strength  via  Weibull  distribution 

The  strength  of  any  material  is  known  to  be  a function  of  somewhat  random 
variables:  manufacturing  flaws,  impurities,  voids,  process  variations,  etc.  Therefore, 
the  strength  of  any  given  material  is  more  accurately  represented  by  a statistical 
distribution  than  by  a single  average  handbook  value.  Freudenthal  [I3J  employed  this 
philosophy  in  his  investigation  of  the  ultimate  strengths  of  aircraft  components 
without  impact  damage.  Various  aircraft  components  and  substructures  were  tested  to 
failure  under  loading  conditions  similar  to  what  they  would  experience  in  operation. 
This  raw  strength  data  was  then  statistically  analyzed  employing  various  probability 
distributions  to  represent  the  observed  scatter.  It  was  found  that  the  Weibull  dis- 
tribution function  [14]  provided  the  best  representation  of  his  residual  strength  data. 
This  distribution  function  is  written  as 

L(z)  = exp  [ -(zM"]  (6) 

where  L(z)  is  the  probability  of  component  survival, 

z = cr/av UT,  the  loading  stress/design  ultimate  stress,  and 
v,  7j  are  the  empirical  parameters  for  Weibull  fit. 

Two  important  facts  are  evident  when  considering  Freudenthal’s  statistical 
results.  First,  the  design  ultimate  stress,  <t-uIT.  is  that  averaged  strength  level  measured 
under  the  considered  type  of  loading.  For  example,  a component  under  tension  would 
be  normalized  to  the  ultimate  tensile  strength  of  the  material,  one  under  torsion -the 
shear  strength,  etc.  Second,  “the  distribution  of  the  ultimate  strength  of  the  test 
specimens  can  be  considered  to  represent  a single  population  irrespective  of  the  type 
of  structures  tested  and  its  mode  of  failure,  as  long  as  this  failure  can  be  classified  as 
ultimate  [13]”.  This  postulate  is  born  of  necessity  to  a degree  because  of  the  cost  of 
replication  and  the  accompanying  scarcity  of  ultimate  strength  data  for  aircraft 
components.  More  data  is  required  to  prove  the  validity  of  the  postulate;  the  resulting 
distribution  obtained  by  Freudenthal  will  be  compared  in  a later  section  to  some 
current  residual  strength  data  for  tension  specimens  containing  no  machined  or  impact 
penetration  damage.  It  provides  a limiting  case  for  the  statistical  distributions  of 
impact  penetration  damaged  components  as  dim  approaches  zero. 

To  evaluate  the  residual  strength  of  impact  penetration  damaged  components,  an 
extension  of  the  Freudenthal  approach  for  undamaged  ones  is  made.  Whereas  in  the 
“undamaged”  case,  all  flaw  sizes  are  of  a microscopic  order  of  magnitude,  i.e.  grain 
sizes,  etc.,  impact  damage  can  introduce  a much  more  extensive  flaw  which  then 
controls  the  residual  strength.  The  resulting  residual  strength  distributions  for  the 
impact  damaged  components  become  a set  of  Weibull  functions;  each  one  established 
for  a given  dim.  The  following  development  outlines  the  statistics  for  determining  the 
Weibull  probability  of  failure  function  from  a given  set  of  residual  strength  data 
collected  from  specimens  with  a constant  dim. 

Given  that  the  probability  of  survival  of  a component  containing  dim  under  an 
operational  loading  stress,  <r,  is  provided  by  (6),  the  probability  of  failure  is 

P(z)=  l-L(z) 

= 1 — exp  [ — (z/*')”]  (7) 
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Note  that  at  this  stage  the  theory  is  strictly  applicable  to  a uniaxial  stress  in  the 
region  of  impact  damage. 

To  isolate  v and  tj,  form  the  expression 

ln  {'"  [i~p(7)]l = ln  (ln  texp  idvy>]} 

= tj  In  (z/y) 

= tj  In  z - rj  In  y (8) 

It  can  be  seen  that  this  last  expression  has  the  form  of  a straight  line,  y = mx  + h, 
where 

x = In  z 

b = — tj  In  v 

m = tj  (9) 

This  relationship  provides  a means  for  the  determination  of  y and  tj  from  a set  of 
k residual  strength  data  specimens.  In  the  current  work  the  data  will  be  the  stresses, 
cr,,  in  uniaxial  tension  at  which  specimens  containing  damage,  dim  - const,  failed  by 
fracture.  From  this  data  the  ratios,  z,  = <tJ<tuvj,  can  be  formed.  The  next  step  is  to 
order  in  ascending  fashion  and  renumber  the  z,’s  such  that  z,  has  the  smallest  and  zk 
the  largest  value. 

For  each  of  the  specimens  the  values  of 

"■'■’MiTd}  (i01 

are  then  calculated,  where  the  cumulative  probability  of  failure  of  any  given  specimen 
at  o-j  is  computed  from  the  finite  data  sample  as 

P(Zi)  = ilk  (11) 

In  order  to  represent  the  data  by  a Weibull  distribution,  the  resulting  test  values, 
(x„  y, ) are  fit  to  a straight  line  using  least  squares.  Note  that  the  definitions  of  z and 
P(z,)  insure  that  the  test  values  will  be  on  a single,  but  perhaps  curved,  line.  The  least 
squares  fit  determines  the  straight  line  which  most  closely  approximates  the  test 
values  line.  For  a set  of  best  fit  parameters  (m0,  b0 ),  Eqns.  (9)  are  used  to  determine 
the  corresponding  Weibull  parameters. 

tj  = m0 

v = exp(-bo/m0)  (12) 

This  algorithm  can  be  easily  programmed  for  a digital  computer.  An  alternate 

means  of  finding  the  best  Weibull  fit  for  a set  of  residual  strength  data  is  through  the 
use  of  probability  paper,  which  in  essence  is  a graphical  procedure  which  duplicates 
the  mathematical  algorithm  presented  above  [15]. 

Close  inspection  of  the  Weibull  distribution  function  leads  to  a physical  in- 
terpretation for  v and  ij.  y can  be  considered  as  indicative  of  the  average  strength 
ratio  for  a given  set  of  residual  strength  data,  while  tj  is  a measure  of  the  dispersion  of 
strength  values.  Large  values  of  tj  indicate  a tight  statistical  distribution  with 
correspondingly  low  scatter. 

Once  the  values  of  y and  tj  have  been  determined  for  a given  material  under 
uniform  tension  and  containing  impact  damage.  dim,  Eqn.  (7)  provides  the  probability 

Ini.  Joum.  of  Fracture.  13  (1977)  409-430 


fc  h 


Probability  based  fracture  mechanic*  for  impact  penetration  damage 


415 


of  failure  for  any  single  value  of  loading  tensile  stress.  The  derivative  function  of  (7) 
is  the  probability  density  function  for  the  Weibull  distribution: 

/(z)  = dP(z)/dz 

= exp  (tj/  i>)(z/  i>  )"''[-  (z/ 1')"]  (13) 

It  is  this  function  which  can  be  used  to  represent  the  strength  distribution  in  a 
conventional  structural  reliability  analysis  where  the  operational  load  is  characterized 
by  a spectrum  rather  than  a single  peak  value  [16].  Figure  3 shows  the  probability  of 
failure  in  this  instance  to  be  the  intersecting  area  between  the  anticipated  operational 
load  distribution  and  the  Weibull  strength  distribution  corresponding  to  the  ap- 
propriate value  of  dim. 


Finally,  in  all  statistical  applications  the  question  arises  as  to  the  degree  which 
the  assumed  distribution  function  fits  a given  set  of  data.  In  addition  to  the  past 
evidence  [13,  14]  which  demonstrates  the  applicability  of  the  Weibull  probability 
function  to  “weakest  link"  type  of  physical  phenomena,  typified  by  impact  damage 
residual  strengths,  a positive  check  should  be  made  for  each  correlation  by  a 
goodness-of-fit  test.  In  the  current  work  the  Chi-squared  test  was  used  to  establish  the 
degree  of  correlation. 

4.  Machined  controlled  damage 

To  gain  insight  into  the  sensitivity  of  the  measures  of  impact  damage  with  respect  to 
the  extreme  forms  of  damage  detail,  sharp  cracks  to  smooth  holes,  and  to  establish  a 
statistical  base  for  residual  strength  under  closely  controlled  damage  for  comparison 
to  the  results  from  tests  with  impact  damage  of  the  same  order,  a series  of  residual 
strength  tests  were  performed  with  7075-T6  aluminum  specimens  as  shown  in  Fig.  2. 
These  tests  were  run  under  uniaxial  tension  and  encompassed  three  specimen 
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thicknesses,  0.050,  0.100,  and  0.250  inch,  to  explore  the  fracture  response  from  near 
plane  stress  to  near  plane  strain  conditions.  The  raw  data  from  these  tests  are 
presented  in  Tables  I to  VI  where  the  data  has  been  grouped  according  to  the  form  of 
damage  detail  and  nominal  values  of  dlm.  Three  primary  forms  of  damage  have  been 
investigated:  (l)  no  damage  where  no  machined  damage  is  present  and  failure  is 
controlled  by  the  inherent  material  weaknesses  which  determine  the  ultimate  tensile 
strength;  (2)  smooth  damage  where  failure  is  initiated  at  a smooth  machined  hole;  (3) 
sharp  damage  where  failure  is  initiated  at  sharp  crack  tips.  Six  nominal  values  of  dim 
were  employed:  0.0";  0.593"  - smooth;  0.160",  0.320",  0.450",  and  0.630"  - sharp. 

All  of  the  cracked  specimens  were  machined  such  that  the  rolling  direction  of  the 
material  was  parallel  to  the  orientation  of  the  crack.  In  this  manner  each  crack  tip  lay 
in  a similar  field  of  anisotropy  with  respect  to  the  material  properties.  The  cracks  were 
produced  by  an  electrical  discharge  machining  process  which  resumed  in  narrow  slots 
approximately  0.006"  wide  and  with  tip  radii  on  the  order  of  0.001".  The  slot  lengths 
were  intentionally  made  0.010"  short  of  the  desired  final  crack  lengths.  The  remaining 
0.0 10" crack  tip  length  was  then  made  under  high  cycle  fatigue  with  the  load  normal  to  the 
crack.  The  fatigue  crack  lengths  were  monitored  by  a closed  circuit  television  system  ; the 
resulting  total  crack  lengths  for  each  specimen  are  reported  as  measured.  A more  detailed 
description  and  assessment  of  the  crack  producing  and  monitoring  process  is  available 
[17]. 

For  both  the  no  damage  specimens  and  those  with  only  smooth  holes,  the  rolling 
direction  was  established  perpendicular  to  the  applied  tensile  load.  The  failure  stress 
in  tension,  <rc,  for  all  specimens  was  determined  from  the  peak  tensile  load  at  fracture 
and  the  actual  width  and  thickness  of  the  specimen,  W and  T. 

The  values  of  crc  from  each  table  were  utilized  along  with  <rULT  as  obtained  from 
the  no  damage  results  of  Table  I for  a statistical  analysis  at  each  nominal  dlm.  Figure  4 
presents  the  probability  of  failure  curves  for  all  of  the  sharp  damage  data  along  with 
the  no  damage  curve  for  comparison.  Note  how  the  values  of  the  Weibull  parameters, 
(p,  tj),  from  the  current  no  damage  results,  (1.005,35.0)  compare  with  Freudenthal's 
[13]  parameters  for  the  undamaged  aircraft  components,  (0.96,  19.0).  Both  dis- 
tributions possess  a v value  which  is  near  unity  since  the  z/s  were  computed  based 
upon  a normalization  by  the  average  undamaged  strengths.  However,  the  tj 


TABLE  1 

Residual  strengths  7075-T6  machined,  controlled  damaged 


Nominal  d,m  = 0.0 

No  damage 

No.  A 

L R W 

T 

0 d,. 

<7,  * 

241 

3.503 

0.048 

0.0 

82073. 

242 

3.503 

0.047 

0.0 

87463. 

243 

3.502 

0.048 

0.0 

85070. 

244 

3.503 

0.047 

0.0 

89285. 

245 

3.502 

0.048 

0.0 

83880. 

246 

3.502 

0.047 

0.0 

83843. 

247 

3.501 

0.048 

0.0 

80929. 

248 

3.503 

0.047 

0.0 

89285. 

249 

3.500 

0.047 

0.0 

85714. 

* »uut  = average  <r,  - 8528? 

t All  dimensions  of  length  in  the  tables  are  inches.  (Refer  to  Fig.  2.) 
All  angles  are  given  in  degrees,  stresses  in  psi.  No  entry  in  a 
column  indicates  no  damage  characterized  by  the  column  heading 
existed  in  the  given  specimen. 
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TABLE  II 

Residual  strengths  707$-T(s  machined,  controlled  damage. 


d.m 

= 0.593 

SMOOTH  DAMAGE 

No. 

A 

L 

R 

W 

T 

/3 

d.m 

crt 

281 

0.2965 

3.504 

0.048 

0.593 

64212. 

282 

0.2965 

3.504 

0.047 

0.593 

65579. 

284 

0.2965 

3.500 

0.048 

0.593 

67857. 

285 

0.2965 

3.500 

0.048 

0.593 

64881. 

286 

0.2965 

3.500 

0.047 

0.593 

67477. 

288 

0.2965 

3.500 

0.048 

0.593 

67262. 

289 

0.2965 

3.500 

0.048 

0.593 

67857. 

209 

0.2965 

3.500 

0.101 

0.593 

68175. 

210 

0.2965 

3.500 

0.101 

0.593 

67893. 

212 

0.2965 

3.500 

0.100 

0.593 

67893. 

213 

0.2965 

3.500 

0.100 

0.593 

69429. 

214 

0.2965 

3.500 

0.100 

0.593 

68571. 

215 

0.2965 

3.500 

0.100 

0.593 

68571. 

217 

0.2965 

3.500 

0.100 

0.593 

68000. 

218 

0.2965 

3.500 

0.100 

0.593 

68286. 

68 

0.2965 

3.505 

0.255 

0.593 

65341. 

69 

0.2965 

3.512 

0.255 

0.593 

64764. 

70 

0.2965 

3.510 

0.255 

0.593 

65024. 

72 

0.2965 

3.505 

0.255 

0.593 

66012. 

73 

0.2965 

3.510 

0.254 

0.593 

64607. 

74 

0.2965 

3.515 

0.255 

0.593 

64262. 

75 

0.2965 

3.510 

0.255 

0.593 

65024. 

76 

0.2965 

3.505 

0.255 

0.593 

67131. 

TABLE  III 

Residual  strengths  7075-T6  machined,  controlled  damage. 


Nominal  d,m  = i 

0.160 

SHARP 

DAMAGE 

No. 

A 

L 

R 

W 

T 

P 

d.m 

(T, 

277 

2.406 

3.500 

0.049 

15. 

0.161 

52478. 

278 

2.375 

3.500 

0.048 

15. 

0.159 

57857. 

297 

0.640 

0.169 

0.1510 

3.500 

0.051 

30. 

0.160 

66106. 

298 

0.620 

0.159 

0.1510 

3.500 

0.050 

30. 

0.155 

65143. 

206 

2.390 

3.500 

0.100 

15. 

0.160 

56571. 

207 

2.375 

3.500 

0.100 

15. 

0.159 

59429. 

226 

0.670 

0.184 

0.15 10 

3.500 

0.100 

30. 

0.168 

65143. 

42 

0.660 

3.500 

0.240 

30. 

0.165 

65000. 

43 

0.650 

3.496 

0.255 

30. 

0.162 

60237. 

65 

2.406 

3.505 

0.254 

15. 

0.161 

51220. 

66 

2.437 

3.495 

0.255 

15. 

0.163 

55542. 

67 

2.437 

3.497 

0.252 

15. 

0.163 

54468. 

97 

0.640 

0.169 

0.1510 

3.501 

0.253 

30. 

0.160 

60965. 

98 

0.655 

0.177 

0.1510 

3.497 

0.254 

30. 

0.164 

63496. 

99 

0.637 

0.168 

0.1510 

3.497 

0.252 

30. 

0.159 

61277. 

143 

2.345 

0.782 

0.3905 

3.498 

0.254 

15. 

0.157 

60102. 

144 

2.366 

0.793 

0.3905 

3.497 

0.251 

15. 

0.158 

62660. 

145 

2.330 

0.775 

0.3905 

3.500 

0.255 

15. 

0.156 

62073. 

147 

2.360 

0790 

0.3905 

3.500 

0.255 

15. 

0.158 

61849. 
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TABLE  IV 

Residual  strengths  7075-7  6 machined,  controlled  damage. 


Nominal  d,m  = 

0.320 

SHARP 

> 

2 

> 

E 

No. 

A 

L 

R 

W 

T 

p 

d.. 

<T, 

274 

1.700 

1.496 

0.05 1 

30. 

0.325 

53170. 

275 

1.265 

3.500 

0.050 

30. 

0.316 

54171. 

205 

1.250 

3.503 

0.101 

30. 

0 712 

48332. 

205 

1.234 

3.504 

0.100 

30. 

0.308 

49087 

33 

0.656 

3.496 

0.241 

45. 

0.328 

52105. 

34 

0.640 

3.505 

0.254 

45. 

0.320 

47851. 

35 

0.656 

.7.495 

0.252 

45 

0.328 

51775. 

36 

0.625 

3.497 

0.253 

45. 

0.312 

50184 

37 

0.655 

3.497 

0.253 

45. 

0.328 

49845 

38 

0.656 

3.495 

0.251 

45. 

0.328 

46851. 

39 

0.679 

3.498 

0.25.3 

45. 

0.340 

48136. 

40 

0.656 

3.500 

0.252 

45. 

0.328 

48299. 

57 

1.255 

3.500 

0.253 

30. 

0.314 

51904. 

58 

1.220 

3.501 

0.253 

30. 

0.305 

48095. 

59 

1.255 

3.503 

0.25 : 

30. 

0.314 

47878. 

60 

1 2.70 

3.500 

0.255 

30. 

0.308 

51092. 

61 

1.235 

3,;oo 

0.255 

30. 

0.309 

49748. 

62 

1.275 

3.509 

0.254 

.30. 

0.319 

46682. 

6.7 

1.290 

3.500 

0.254 

.70. 

0.322 

49494. 

64 

1.240 

3.496 

0.255 

30. 

0.310 

50590. 

94 

0.712 

0.145 

0.211 

3.503 

0.251 

45. 

0.356 

50497. 

95 

0.686 

0.132 

0.211 

3.500 

0.253 

45. 

0.343 

54884. 

96 

0.680 

0.129 

0.211 

3.501 

0.252 

45. 

0.340 

52139. 

128 

1.320 

0.262 

0.398 

3.501 

0.253 

30. 

0.330 

56901. 

129 

1 228 

0.216 

0.398 

3.500 

0.249 

30. 

0.307 

61274. 

151 

1 .230 

0.217 

0.398 

3.500 

0.250 

30. 

0.308 

62857. 

132 

1,225 

0.215 

0.398 

3.501 

0.251 

30. 

0.306 

60313. 

133 

1.242 

0.098 

0.523 

3.500 

0.253 

30. 

0.310 

56917. 

134 

1.227 

0.090 

0.52.7 

3.500 

0.25.7 

.70. 

0.307 

56917. 

136 

1.227 

0.090 

0.523 

3.500 

0.254 

30. 

0.307 

57368. 

137 

1.244 

0.099 

0.523 

3.500 

0.252 

30. 

0.311 

56463. 

parameters  indicate  a much  larger  scatter  in  the  Freudenthal  test  data  than  was 
present  in  the  controlled  no-damage  tensile  tests.  The  larger  scatter  in  Freudenthal’s 
data  may  be  attributed  to  his  postulate  for  combining  residual  strengths  from  failure 
tests  under  different  kinds  of  loading. 

One  can  see  from  Fig.  4 that  a structural  engineer  can  assess  the  probability  of 
failure  of  a component  with  machined  damage  given  an  estimate  of  the  principal 
loading  stress  and  the  designed-for  damage  measure.  d:„.  Note  that  as  the  degree  of 
damage  increases  the  probability  curves  shift  towards  the  left.  This  is  consistent  with 
the  concept  that,  at  a given  load  level,  the  probability  of  failure  increases  with  the 
degree  of  damage;  a concept  which  will  be  seen  to  carry  over  into  impact  damage 
probability  curves. 

Figure  5 shows  how  the  type  of  damage  detail  is  reflected  through  the  probability 
curves.  The  two  groups  have  roughly  the  same  dim  where  in  fact  the  smooth  damage  is 
greater  than  the  sharp.  However,  the  sharp  crack  controlled  damage  curve  is  located 
significantly  to  the  left  of  the  smooth  stress  concentration  controlled  damage  curve. 
This  demonstrates  the  importance  of  one  aspect  of  damage  detail,  smoothness  versus 
sharpness,  and  is  useful  in  interpreting  the  scatter  of  actual  impact  damage  residual 
strengths.  The  degree  to  which  any  given  set  of  impact  damage  patterns  is  charac- 
terized by  smoothness  or  sharpness  depends  upon  the  exact  damage  detail  and  is 
therefore  directly  influenced  by  the  impact  conditions  including  the  energy  level  of  the 

I impact  phenomenon. 

The  effect  of  specimen  thickness  is  also  shown  by  a shift  in  the  probability 

i I 
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TABLE  V 

Residual  strengths  7075-76  machined . controlled  damage. 


Nominal  d,m  = 

0.450 

SHARP 

DAMAfi 

E 

No. 

A 

L 

R 

W 

T 

P 

d„ 

(T, 

265 

0.875 

3.499 

0.051 

45. 

0.438 

53797. 

267 

0.880 

3.498 

0.050 

45. 

0.440 

55575. 

268 

0.880 

3.497 

0.050 

45. 

0.440 

55591. 

269 

0.890 

3.497 

0.049 

45. 

0.445 

56025. 

270 

0.895 

3.500 

0.050 

45. 

0.448 

54514 

271 

0.885 

3.4% 

0.049 

45. 

0.442 

56041 

272 

0 895 

3.4% 

0.050 

45. 

0.448 

52860. 

273 

0.880 

3.4% 

0.050 

45. 

0.440 

54062. 

293 

0.640 

0 086 

0.234 

3.500 

0.050 

60. 

0.480 

49143 

293 

0.640 

0.086 

0.234 

3.500 

0.049 

60. 

0.480 

52478. 

194 

0.895 

3.497 

0.099 

45. 

0.448 

44194. 

195 

0.890 

3.500 

0.100 

45. 

0.445 

43714. 

1% 

0.885 

3.500 

0.100 

45. 

0.442 

43429. 

197 

0.875 

3.497 

0,099 

45. 

0.438 

44194. 

198 

0.885 

3.4% 

0.100 

45. 

0.442 

45767. 

199 

0.875 

3.497 

0.101 

45. 

0.438 

44168. 

200 

0.890 

3.500 

0. 101 

45. 

0.445 

43281. 

201 

0.885 

3.500 

0.099 

45. 

0.442 

45310. 

222 

0.635 

0.084 

0.234 

3.500 

0.099 

60. 

0.476 

45.310. 

224 

0.625 

0.079 

0.234 

3.500 

0.099 

60. 

0.469 

47619. 

31 

0.656 

3.4% 

0.254 

60. 

0.492 

40879. 

54 

0.968 

3.505 

0.254 

45. 

0.484 

40437. 

55 

0.937 

3.497 

0.253 

45. 

0.468 

40916. 

56 

0.930 

3.500 

0.253 

45. 

0.465 

39977. 

123 

0.890 

0.211 

0.234 

3.500 

0.252 

45. 

0.445 

44218. 

124 

0.880 

0.206 

0.234 

3.502 

0.252 

45. 

0.440 

44872. 

125 

0.880 

0.206 

0.234 

3.500 

0.253 

45. 

0.440 

44043. 

126 

0.906 

0.219 

0.234 

3.502 

0.253 

45. 

0.453 

41309. 

curves.  Figure  6 illustrates  the  probability  curves  computed  for  the  data  of  Table  VI 
when  grouped  according  to  specimen  thickness.  The  shift  of  the  probability  curves  to 
the  right  for  the  thinner  specimens  could  be  ascribed  to  a plane  stress  effect  on 
residual  strengths.  However,  the  relatively  low  number  of  specimens  at  the  thinner 
thicknesses  and  the  overlapping  of  the  curves  at  low  probability  levels  precludes  any 
strong  conclusions  in  support  of  a large  thickness  effect  on  average  strengths.  For  the 
range  of  thicknesses  investigated,  data  from  all  thicknesses  at  a given  nominal  du* 
were  grouped  together  for  the  purpose  of  generating  the  probability  curves  of  Fig.  4. 

The  manner  of  evaluation  of  the  goodness-of-fit  of  the  Weibull  distribution  to 
these  data  groupings  via  the  Chi-squared  test  is  described  in  the  appendix.  The 
correlation  coefficient  for  each  fit  is  given  with  the  probability  curves.  Note  that  a 1.00 
correlation  indicates  a perfect  Weibull  fit  via  a Chi-squared  measure.  Interpretation  of 
the  relatively  low  correlation  numbers  obtained  for  the  current  data  is  discussed  in  a later 
section. 


5.  Impact  damaged  specimens 

As  an  example  of  the  application  of  probability  based  fracture  mechanics  for  impact 
penetration  damage,  the  following  statistical  analysis  was  performed  with  residual 
strength  data  reported  in  the  literature  [8,  18,  19].  In  these  references  residual  strength 
tests  were  conducted  upon  high  strength  aluminum  panels  that  were  first  subjected  to 
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TABLE  VI 

Residual  strengths  707 V7T>  machined,  controlled  damage 


Nominal  </,„  = 

0.630 

SHARP 

DAMAGE 

No 

A 

L 

R 

W 

T 

0 

d. 

(T, 

250 

0.630 

3.500 

0.047 

90. 

0.6.30 

43161 

251 

0.620 

3.500 

0.047 

90. 

0.620 

42310. 

252 

0.620 

3.500 

0.047 

90. 

0.620 

42310. 

253 

0.625 

3.500 

0.047 

90. 

0.625 

42553. 

254 

0.615 

3.500 

0.047 

90. 

0.615 

43465. 

255 

0.615 

3.500 

0.047 

90. 

0.615 

410.33. 

256 

0.620 

3.500 

0.047 

90. 

0.620 

38298. 

259 

0.656 

3.495 

0.050 

75. 

0.612 

53791. 

260 

0.650 

3.500 

0.05 1 

75. 

0.606 

53109. 

290 

0.642 

0.122 

0.199 

3.500 

0.048 

90. 

0.642 

38571. 

181 

0.615 

3.502 

0.099 

90. 

0.615 

41823. 

182 

0.625 

3.501 

0.099 

90. 

0.625 

40392. 

183 

0.635 

3.500 

0.099 

90. 

0.635 

41847. 

185 

0.610 

3.501 

0.099 

90. 

0.610 
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Figure  4 Probability  curves  for  high  strength  aluminum,  machined,  controlled  damage  specimens 
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Figure  5.  Influence  of  damage  type  on  residual  strength  probability  curves. 
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Figure  6.  Thickness  effect  on  residual  strength  probability  curves  for  machined,  controlled  damage 
specimens. 

ballistic  impact  damage  and  then  pulled  to  fracture  in  a uniaxial  tensile  test.  The  raw 
data  from  these  tests  were  sorted  and  assembled  into  Table  VII.  Note  that  in  this  table 
the  results  were  ordered  by  dm.  The  data  were  originally  obtained  in  a manner 
designed  to  study  a spectrum  of  impact  damage  sizes  over  a range  of  materials  and 
specimen  thicknesses.  For  correlation  with  the  Weibull  function  and  comparison  to 
the  results  of  the  machined,  controlled  damage  specimens,  Table  VII  is  comprised 
only  of  data  which  was  taken  on  aluminum  specimens  ai  thicknesses  near  those  of  the 
controlled  damage  tests.  Also  residual  strength  data  was  chosen  from  the  references 
only  in  cases  where  dlm  could  be  reliably  computed. 

Because  of  the  spread  in  d ^ values,  the  data  was  grouped  for  two  nominal 
values,  1.70  and  2.60  inches,  within  ±20%  in  each  case.  The  residual  strength  values 
were  normalized  by  an  ultimate  strength  value  consistent  with  the  reported  material 
properties. 

Figure  7 presents  the  resulting  probability  of  failure  curves  for  impact  damage  in 
high  strength  aluminum. 


6.  Discussion  of  results 


Figure  8 presents  v and  i\  versus  d ^ for  each  set  of  residual  strength  data  analyzed.  It 
can  serve  to  focus  several  key  concepts  that  have  evolved  during  this  current  work. 
First,  again  note  that  at  no  impact  damage,  dm  = 0.0,  the  values  for  v and  t)  are  given 
for  the  machined  aluminum  specimens  and  the  Freudenthal  work  on  aircraft  structural 
components.  Then  for  the  controlled  machined  damaged  cracked  specimens,  the 
average  strength  values  as  reflected  by  v fall  as  d*  is  increased.  The  corresponding 
scatter  of  residual  strengths  is  seen  to  rise  somewhat  as  indicated  by  a drop  in  tj  over 
the  same  range  of  dm.  It  is  interesting  to  observe  the  measurable  increase  in  v for  the 
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TABLE  VII 

Residual  strengths  impact  damage. 


A:  7075-T6  Aluminum,  rolling  grain  transverse  to  loading  axis 

B:  7075-T6  Aluminum,  rolling  grain  parallel  to  loading  axis 

C:  2024-T8I  Aluminum,  rolling  grain  transverse  to  loading  axis 

D:  2024-T3  Aluminum,  rolling  grain  transverse  to  loading  axis 

Materia) 

T 

W 

P* 

0 

d„ 

<Jr  0 

'J  O'lI  t 

REF. 

(in) 

(in) 

(in) 

O 

(in) 

(psi) 

Nominal  d, 

- « 

1.70  (± 

20%) 

C 

0.160 

12. 

1.35 

80. 

1.31 

41800. 

0.624 

(81 

A 

0.032 

12. 

1.45 

80. 

1.41 

41000. 

0.539 

[81 

B 

0.250 

12. 

1.45 

90. 

1.45 

60700. 

0.720 

[191 

B 

0.500 

12. 

1.47 

90. 

1.47 

56500. 

0.670 

[19] 

A 

0.032 

12. 

1.53 

90. 

1.53 

37800. 

0.497 

[8| 

A 

0.250 

12. 

1.58 

90. 

1.58 

49400. 

0.648 

18) 

B 

0.500 

12. 

1.58 

90. 

1.58 
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(19) 

A 

0.500 

7. 

1.65 

80. 

1.60 

45000. 

0.592 

[8] 

B 
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12. 

1.61 

90. 

1.61 
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0.590 

119) 

B 

0.125 

12. 

1.65 

85. 

1.64 

58600. 

0.773 

(8) 

B 

0.500 

8. 

1.70 

90. 

1.70 
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0.720 

[191 

B 

0.190 

18. 

1.70 

90. 
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0.546 
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A 
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70. 
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A 
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7. 
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0.610 
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B 
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0.630 

[19] 

D 
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90. 
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0.692 
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A 

0.250 

12. 

1.90 
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B 
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1 92 
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1.95 
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1.95 
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0.490 
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» = 
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B 
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12. 

2.14 

90. 

2.14 
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[19] 

A 
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12. 

2.16 

90. 

2.16 

41200. 

0.542 

[8] 

A 

0.125 

12. 

2.30 

90. 
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34000. 

0.446 

[8] 

B 

0.125 

12. 

2.35 

90. 
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43300. 

0.570 

[8] 

A 
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9. 

2.40 

90. 
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52000. 

0.686 

[8] 

B 

0.250 

12. 

2.42 

90. 

2.42 

45400. 

0.540 

[19] 

A 

0.190 

18. 

2.96 

64. 

2.39 

30700. 

0.405 

[18] 

A 

0.190 

18. 

2.96 
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2.47 

21500. 

0.284 

[18] 

B 
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12. 

2.56 

90. 

2.56 
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0.630 

[19] 

A 

0.090 

12. 

2.70 

90. 

2.70 

29400. 

0.387 
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C 

0.250 

12. 

2.75 

85. 

2.73 

27500. 

0.410 

[8] 

B 

0.375 

12. 

2.99 

75. 

2.79 
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0.360 

119) 

B 
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16. 

2.82 

90. 

2.82 

40100. 

0.470 
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B 

0.250 

16. 

2.85 
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2.85 
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0.570 

119] 

C 

0.250 

12. 

3.22 

80. 

3.12 

23400. 

0.350 

[8] 

C 

0.250 

12. 

3.25 

80. 

3.15 

20400. 

0.305 

[8] 

* When  photograph  of  damage  was  not  given  in  reference,  tabulated  LD  for  high 
angles  of  0 was  assumed  to  he  reliable  estimation  of  p. 
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Figure  7.  Probability  curves  for  high  strength  aluminum,  impact-damaged  specimens. 
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Figure  8.  Weibull  parameters  vs.  impact  damage  measure. 
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machined  hole  specimens  over  the  cracked  ones  along  with  an  tj  value,  which  was 
some  five  to  ten  times  that  for  the  cracked  damage.  Finally,  note  the  apparent  increase 
in  average  strengths  for  the  impact  damaged  specimens  over  the  trend  in  v established 
by  the  controlled  damage  specimens.  This  was  accompanied  by  a continued  increase 
in  scatter  as  indicated  by  the  further  decline  in  tj  values.  After  consideration  of 
photographs  of  the  impact  damaged  specimens  reported  in  Ref.  [19],  it  was  concluded 
that  both  the  increase  in  scatter  and  apparent  rise  in  the  average  strength  of  the 
impact  damage  results  were  due  in  part  to  the  nature  of  the  damage  detail.  For 
example,  for  many  of  the  thicker  specimens,  the  high  velocity  projectiles  used  to 
penetrate  the  panels  often  created  damage  which  consisted  primarily  of  a through  hole 
plus  an  extended  spall  surface  on  the  exit  side.  Little  evidence  of  extensive  radial 
crack  formation  was  indicated.  Therefore,  those  specimens  acted  more  like  the 
smooth  damage  of  the  controlled  damage  tests  and  caused  an  apparent  increase  in  the 
average  residual  strength.  This  fact  may  also  account  for  an  apparent  reversal  of  the 
plane  stress  toughness  effect  for  impact  residual  strengths  as  reported  in  the  literature 
[20]. 

The  residual  strengths  for  the  smooth  damage  specimens  had  the  least  scatter  as 
measured  by  tj  of  any  of  the  test  groups.  It  was  better  than  the  standard  tensile  test 
results  by  approximately  30%.  This  implies  that  the  stress  concentration  in  the  smooth 
hole  specimens  strongly  focused  the  conditions  for  failure  such  that  the  variation  in 
resulting  residual  strengths  was  less  than  that  inherent  for  failure  under  gross  necking 
as  associated  with  the  uniaxial  tension  test.  Also,  the  specimen  thickness  appears  to 
have  had  a minimal  influence  in  the  scatter  of  the  smooth  damage  data.  In  contrast  the 
variability  of  the  sharp  machined  damage  residual  strengths  is  measurably  higher  as 
indicated  by  their  lower  tj  values.  This  could  reflect  the  sensitive  nature  of  the 
mechanism  of  crack  tip  extension  to  local  crack  tip  geometry  and  specimen  thickness. 
In  fact.  Fig.  6 shows  that  in  the  case  of  the  constant  thickness,  0.250  inch,  specimens 
with  crack  damage,  d = 0.630,  a Weibull  analysis  on  thirty-three  specimens  resulted  in 
a significantly  reduced  scatter  in  residual  strengths  over  the  entire  group  of  Table  VI 
which  included  all  three  thicknesses  at  a dim  = 0.630;  tj  = 29.03  versus  10.7.  These 
points  become  important  when  interpreting  impact  damage  results  where  a mixture  of 
smooth  and  sharp  damage  detail  produces  higher  average  strengths,  similar  to  smooth 
damage  data,  but  significantly  lower  tj  values,  typical  to  the  variability  influences 
associated  with  sharp  damage  data. 

Whereas  v and  tj  reflect  the  average  strength  and  scatter  respectively  for  a given 
data  set,  the  Chi-squared  correlation  number  reflects  the  degree  to  which  the  Weibull 
density  function's  shape  corresponds  to  that  of  the  histogram  of  experimental  cu- 
mulative probabilities.  The  following  appendix  reveals  in  detail  that  for  the  data  under 
current  consideration,  two  separate  cases  arise  with  respect  to  the  low  correlation 
numbers.  For  the  case  of  the  machined,  controlled  damage  specimens,  the  low 
correlation  number  is  due  to  the  non-random  nature  of  each  test  data  sample  as  a 
result  of  the  sub  set  groupings  at  different  but  specified  thicknesses.  However,  the  low 
correlation  number  has  little  practical  significance  because  the  resulting  probability 
curves  of  Figs.  4 to  6 are  each  the  integral  of  a density  function  which  has  its 
maximum  influence  over  a small  interval  of  the  z domain.  It  is  therefore  insensitive  to 
the  exact  shape  of  the  density  function  and  depends  more  upon  the  average  strength 
and  scatter  parameters. 

For  the  case  of  impact  damage,  the  statistical  scatter  extends  over  a much  greater 
portion  of  the  z domain,  and  the  correlation  number  becomes  more  important  in 
reflecting  possible  significant  error  in  the  resulting  probability  curves.  The  correlation 
numbers  of  Fig.  7 for  the  impact  damage  probability  curves  are  seen  to  be  sig- 
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nificantly  improved  over  the  controlled  damage  case.  In  fact,  given  that  the  number  of 
impact  damage  specimens  available  for  statistical  analysis  was  near  the  minimum 
allowable  for  a meaningful  Chi-squared  test,  the  resulting  correlation  numbers  are 
perhaps  as  good  as  can  be  expected. 

These  results  point  to  the  need  for  careful  consideration  of  the  following 
guidelines  for  the  future  application  of  a probability  based  fracture  mechanics 
approach  to  impact  penetration  damage. 

i.  Isolate  and  group  specimens  according  to  material,  thickness,  and  dim  in  order 
to  derive  the  most  reliable  statistics  as  measured  by  a maximum  possible  tj  value  and 
a maximum  Chi-squared  correlation. 

ii.  Test  with  enough  replications  at  each  grouping  to  ensure  a meaningful  cor- 
relation via  a Chi-squared  test.  A minimum  number  of  replications  appears  to  be 
twenty  to  thirty;  a better  number  would  be  near  fifty. 

Because  of  the  cost  factor  associated  with  these  recommendations,  future 
research  into  the  development  of  simulation  techniques  via  numerical  computer  to 
obtain  probability  of  failure  curves  for  impact  damage  would  appear  to  offer  a 
potential  alternative  to  minimize  expense. 

7.  Conclusions 

This  work  represents  the  first  step  in  the  development  of  a statistical  approach  for 
evaluating  the  residual  strength  of  impact  penetration  damaged  structures.  It  es- 
tablishes a philosophy  and  accompanying  methodology  for  the  assessment  of  the 
degree  of  damage  resulting  from  an  impact  penetration  along  with  corresponding 
evaluation  of  the  residual  strength  in  terms  of  a probability  of  failure.  This  is  in 
contrast  to  previous  fail/no  fail  deterministic  approaches  which  are  unable  to  account 
for  large  scatter  in  the  residual  strength  data  taken  from  impact-damaged  specimens. 
As  such,  the  statistical  approach  is  of  immediate  use  to  the  engineering  community  in 
beginning  to  include  quantitative  probability  based  fracture  mechanics  in  early  phases 
of  design  for  impact  penetration  damage.  As  a larger  data  base  is  collected  following 
the  guidelines  recommended  in  this  paper,  precise  calculations  of  the  probability  of 
failure  for  impact-damaged  structures  will  be  possible  beyond  the  estimations  now 
attainable  for  the  limited  damage  range  of  Figs.  7 and  8. 

Recognizing  that  the  formulation  and  parameters  presented  in  this  paper  are 
strictly  for  the  case  of  uniaxial  stress  states,  future  research  will  be  involved  in 
generalizing  the  statistical  approach  for  multiaxial  stress  fields.  Note  that  past  workers 
[81  have  attempted  to  extend  uniaxial  test  data  by  characterizing  the  stress  field 
remote  from  the  damage  in  terms  of  the  largest  principal  stress.  While  such  an 
approach  could  be  used  with  the  probability  functions  as  provided  in  this  paper,  it  is 
recommended  only  in  cases  where  one  stress  component  obviously  dominates  the 
stress  field  in  the  vicinity  of  the  damage. 

For  more  general  states  of  multiaxial  stress,  two  avenues  of  research  must  be 
explored.  One  consists  of  seeking  ways  to  combine  the  uniaxial  probabilities  to  form  a 
practical  estimation  of  the  probability  of  failure  for  an  impact  penetration  damage  in  a 
multiaxial  stress  field.  The  second  involves  the  development  of  a multiaxial  version  of 
dim  with  associated  probability  correlations.  For  example,  following  the  manner  of 
derivation  of  d,m  given  in  this  paper,  but  where  the  inclined  crack  is  situated  in  a 
multiaxial  stress  field.  <r,  and  <r2.  a relationship  between  the  fracture  stresses  and  the 
damage  parameters,  p and  /3.  becomes 

C«=  <rf(irp/2)[sin2  /3  + (o-j/o-,)2  cos2  p + 4(<r2/<r,)  sin2  /3  cos2  f$\  (14) 

where  /3  is  measured  from  the  direction  of  a-,. 
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This  expression  illustrates  that  the  damage  factor  is  not  a simple  function  of  p 
and  p as  in  the  uniaxial  case,  but  contains  the  ratio  of  the  multiaxial  stress 
components  as  well.  More  work  is  obviously  required  to  evaluate  this  and  alternative 
precise  measures  of  impact  damage  in  a multiaxial  stress  field. 

Finally,  recall  that  the  Weibull  two-parameter  probability  function  was  utilized  in 
this  current  effort  based  upon  the  earlier  work  by  Freudenthal,  who  identified  it  as  the 
best  function  to  use  in  the  statistical  representative  of  residual  strengths  of  un- 
damaged structural  components.  It  is  important  to  note,  however,  that  the  probability 
based  fracture  mechanics  approach  is  not  limited  to  any  one  density  function.  Future 
work  should  explore  the  utilization  of  other  density  functions  to  represent  the  scatter 
of  residual  strengths  of  impact-damaged  components  with  comparisons  based  upon 
the  four  factors  for  goodness  of  fit  assessment  outlined  in  the  appendix  of  this  paper. 

Appendix 

The  Chi-squared  test  for  goodness  of  fit  provides  a means  for  estimating  how  closely 
a given  Weibull  function  matches  the  distribution  of  residual  strength  data  from  which  it 
was  derived.  Figure  9 illustrates  the  distribution  of  residual  strengths  in  the  form  of  a 
histogram  for  the  test  data  of  Tables  II  and  III  where  d,m  =0.593  and  0.160  inch, 
respectively.  Superimposed  upon  each  figure  is  the  corresponding  Weibull  probability 
density  function.  The  ordinate  scale  has  been  chosen  such  that  the  total  area  under  the 
density  function  is  unity  and  the  histogram  was  constructed  such  that  the  area  which  it 
encloses  is  also  equal  to  unity. 

A key  step  in  establishing  a correlation  number  by  the  Chi-squared  test  is  the 
division  of  the  probability  density  area  into  known  probabilities;  it  is  convenient  to 
make  these  sub-areas  equal.  The  intervals  on  the  z axis  which  establish  the  boun- 
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Figure  9.  Chi-squared  correlation  for  Weibull  fits. 
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daries  of  the  sub-areas  are  determined  from  the  inversion  of  the  Eqn.  (7). 

z,=  ir[  — ln(l  - (Al) 

where  P(z()  is  the  cumulative  probability  up  to  and  including  interval,  j. 

Comparison  is  then  made  between  the  number  of  observed  test  values  which  lie 
in  each  sub-area  and  those  that  would  be  expected  based  upon  the  sub-area  prob- 
ability and  the  total  number  of  test  strength  values  in  the  sample.  This  comparison  is 
made  through  the  Chi-squared  parameter  defined  by 

N 

X1  = 2 (O,  - EflE;  (A2) 

/-i 

where 

E,  = expected  number  in  interval,  j 

O,  = observed  number  in  interval,  j 

N = total  number  of  sub-area  intervals 

Based  upon  the  Chi-squared  value  and  the  number  of  independent  sub-area 
intervals,  tables  are  available  [21]  which  indicate  the  degree  of  correlation  between 
the  Weibull  distribution  and  the  actual  residual  strength  distribution.  The  number  of 
independent  sub-areas  are  termed  the  degrees  of  freedom  and  for  the  two-parameter 
Weibull  fit  is  given  by 

dof  = N-  3 (A3) 

where  the  constant,  3.  accounts  for  the  dependent  relationship  between  the  ex- 
perimental data  and  the  two  Weibull  parameters  along  with  the  constraint  that  all  the 
sub-area  probabilities  must  add  to  unity. 

From  consideration  of  Fig.  9,  four  factors  emerge  as  equally  important  in 
assessing  how  well  a given  Weibull  fit  corresponds  to  a distribution  of  residual 
strengths  and  how  well  a given  Weibull  density  function  may  represent  the  actual 
probability  of  failure  curve.  First,  the  shape  and  relative  size  of  the  histogram  and 
density  function,  when  both  are  normalized  to  the  same  total  area,  provide  a check  of 
the  suitability  of  the  assumed  distribution.  For  the  case  of  the  0.160  data  the  Weibull 
distribution  appears  to  correspond  quite  well  to  the  normalized  histogram.  However,  a 
detailed  plot  of  the  histogram  and  Weibull  function  for  the  0.593  data  shows  that  the 
actual  strength  data  reflect  a bi-modal  distribution  which  can  be  traced  to  a thickness 
effect  in  the  residual  strengths  of  panels  with  smooth  holes  and  which  is  not 
accounted  for  by  a Weibull  distribution. 

Second,  the  Chi-squared  value  itself  is  greatly  influenced  by  the  total  number  of 
specimens.  For  example  in  both  cases  shown  on  Fig.  9.  the  total  number  of  specimens 
in  each  set  was  such  that  the  expected  value  in  each  interval  was  some  fractional 
number  near  five.  Thus,  since  the  observed  value  in  each  interval  must  be  an  integer 
near  five,  an  error  arises  from  the  fractions  which  can  significantly  alter  *2.  In  the 
case  of  the  0.160  data,  this  along  with  the  larger  difference  in  interval  2 has  lead  to  the 
rather  low  correlation  number,  0.142.  While  this  third  factor,  the  correlation  number, 
itself,  is  an  important  number  in  strictly  defining  a correspondence  between  statistical 
fit  and  experimental  data,  it  is  highly  sensitive  to  inaccuracies  arising  from  a low 
number  of  specimens.  For  practical  purposes,  the  density  function  having  a low 
correlation  number  may  yet  yield  a very  useful  probability  of  failure  curve. 

The  0.593  data  is  a case  in  point  as  shown  in  Fig.  9 and  illustrates  the  fourth 
factor  in  assessing  statistical  fits.  Here  it  is  obvious  that  the  detailed  bi-modal  nature 
of  the  data  cannot  be  closely  fit  by  the  Weibull  distribution,  and  a very  low 
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correlation  number  results.  However,  because  the  distribution  is  very  narrow,  the 
bi-modal  nature  of  the  sample  is  a higher-order  effect  with  respect  to  the  resulting 
shape  and  location  of  the  probability  of  failure  curve.  The  Weibull  function  for  that 
data  set  accurately  reflects  both  the  average  strength  and  scatter  band.  It  leads  to  a 
probability  of  failure  curve  which  is  quite  accurate  and  useful  for  practical  purposes 
even  though  the  correlation  number  is  low.  Therefore,  assessment  of  the  effect  of  the 
width  of  the  data’s  scatter  band  upon  the  computational  accuracy  of  the  probability 
curve  is  an  important  factor  in  evaluating  the  need  for  a close,  detailed  matching  of 
the  experimental  histogram  and  the  derived  Weibull  density  function. 

In  conclusion,  all  four  factors  become  important  in  assessing  the  degree  of 
success  in  representing  the  statistical  nature  of  a given  set  of  test  data.  In  order  to 
optimize  the  goodness-of-fit  by  all  four  factors  in  future  extensions  of  the  statistical 
approach,  the  guidelines  presented  in  the  discussion  section  of  this  paper  are  recom- 
mended. 
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RESUME 

Le  memoire  presente  une  nouvelle  methode  pour  verifier  la  resistance  residuelle  de  composants  structuraux 
ayant  subi  un  dommage  suite  4 un  impact  penetrant.  La  methode  est  basee  sur  ('application  de  la  fonction  de 
density  de  probability  de  Weibull,  pour  tenir  compte  de  la  dispersion  appreciable  de  resistances  residuelles,  qui 
decoule  de  la  nature  statistique  du  dommage  par  impact  (fissure,  trous  et  ecaillages  en  surface). 

Pour  mesurer  l'etendue  d'une  configuration  d'un  dommage  par  impact,  on  presente  un  nouveau  parametre 
base  sur  la  solution  applicable  a une  fissure  inclinee  dans  un  champ  de  tension  uniaxial,  et  base  sur  un  critere  de 
rupture  obtenu  par  voie  experimentale,  pour  des  configurations  de  fissure  de  mode  mixte. 

Les  resultats  des  essais  sont  presentes  pour  des  resistances  residuelles  d'eprouvettes  endommagees  et 
usinees  afin  de  verifier  la  nature  statistique  qui  caracterise  un  detail  d endommagement  exterieur,  un  trou  a 
paroi  lisse  et  une  fissure  aigue.  Enfin,  des  courbes  de  probability  de  rupture  sont  fournies  pour  des  panneaux 
d'aluminium  7075-TG  soumis  a tension  uniaxiale  et  comportant  a la  fois  un  endommagement  par  usinage  et  par 
impact. 
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